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Abstract: In this paper we present a two channel model with the goal of reproducing the soft scattering 
data available in the ISR-Tevatron energy range, and extend the model results to LHC and Cosmic Rays 
energies. A characteristic feature of the model is that we represent the sum of all diffr active final states at 
a vertex, by a single diffractive state. Our two main results are: (i) The approach of the elastic scattering 

' J amplitude to the black disc bound is very slow, reaching it at energies far higher than the GZK ankle 
cutoff, (ii) Our predicted survival probability for Higgs exclusive central diffractive production at the LHC 
is 0.7%, which is considerably smaller than our previous estimate. The above features are compatible 

•O with a parton-like model in which the traditional soft Pomeron is replaced by an amplitude describing the 
h> partonic system, which is saturated in the soft (long distance) limit. 
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1. Introduction 

The GLM model is based on a two channel eikonal like solution of the s-channel unitarity equation [1]. 
Our present investigation is based on an improved model, and we calculate and discuss it's predictions 
and implications at the LHC and Cosmic Rays energies. For the correct degrees of freedom, the partial 
amplitude for scattering of state %' with state '£;' can be written as 

4f(s,b) = i8 i>v 8 k , k , (l - exp . (1.1) 
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i and k represent sets of quantum numbers which eliagonalize the interaction matrix. The opacities fij^ 
are arbitrary real functions of energy and impact parameter. 

As long as are not explicitly specified, our presentation is model independent. We assume that the 
scattering amplitude at high energies is predominantly imaginary. In our original single channel model [2] 
we assumed that the opacity is determined by the exchange of a soft Pomeron represented by a factorized 
fixed pole in the complex angular momentum (J) plane. This simplified assumption is not maintained in 
the present, more elaborate two channel model, where factorization of the coupling constants is relaxed. 
These features are compatible with our present partonic picture, where the soft Pomeron, is replaced by 
the soft distance limit, of the amplitude for a saturated partonic system [3]. This approach is supported 
by eikonal-like models [4,5], which reproduce the e-p DIS data over a wide range of Q 2 , starting from very 
small virtualities. Details of our parametrization are presented in Sec. 2. 

We shall elaborate on the following issues: 

1) Our original investigation [1] neglected the double diffraction state. Experimental data on this channel 
is now available [6] and it makes the simplified two amplitude approximation doubtful. The present three 
amplitude analysis is based on an updated data base which includes the published p-p and p-p data points 
of (Jtoti the integrated values of a e i, a s d, (Jdd an d the forward elastic slope B e \ in the ISR-Tevatron energy 
range. The forward slopes of the SD and DD final states, as well as p = /^""'/(t^o's) ' are predictions of the 
model. 

2) Based on this initial investigation, our present analysis aims at providing reliable predictions of the 
quantities to be measured at 14 TeV, the LHC cm. energy. Our output also covers the broad Cosmic 
Ray energy range up to the GZK limit. We validate the consistency of our calculations with unitarity by 
extending our output up to the Planck mass. 

3) One requires a reliable formulation of soft scattering to calculate the survival probability of large rapidity 
gaps (LRG), initiated by the underlying soft rescatterings of the spectator partons, in an inelastic diffractive 
(soft or hard) process [2,7,8]. This calculation is of particular importance for the assessment of the 
discovery potential for LHC Higgs production in an exclusive central diffractive process. This channel, 
with a clean two LRG signature, has a relatively good signal to background ratio. The extraction of a 
clear diffractive Higgs signal at the LHC, requires a precise knowledge of the cross section and transverse 
momentum behavior of the single diffractive production channel, which provides a significant contribution 
to the background of interest. We shall investigate the SD channel in some detail. 

4) Some of the fundamental consequences of s-channel unitarity in the high energy limit are not clear 
as yet. We wish to assess the rate that the elastic scattering amplitude reaches the unitarity black disc 
bound in 6-space. This is, obviously, coupled to the behavior of the corresponding diffractive amplitudes. A 
related piece of information, which is still unknown, is the rate at which the proton black core expands with 
energy in 6-space. At present, different models provide drastically different assessment of this phenomenon. 
Our advantage of having a specific model enables us to study in detail the behavior of the amplitudes as 
functions of energy and impact parameter, and, consequently, we are able to provide a numerical description 
of the process. Our approach differs from alternative treatments which are based on general assessments 
which are not accompanied by a detailed analysis of the soft scattering data. See, for example Ref. [9]. 
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The organization of this paper is as follows: In Sec. 2 we briefly summarize the general properties 
of the eikonal approach and formulate our model. Following, we review the extension from GLM single 
channel to two channels. In Sec. 3 we present the three models we have considered and the output of our 
calculations in the ISR-Tevatron energy range. Sec. 4 is devoted to a presentation of our predictions for 
LHC and Cosmic Ray energies. In Sec. 5 we discuss survival probability calculations in the GLM models 
applied to Higgs production at LHC. Sec. 6 is devoted to a detailed analysis of the onset of unitarity effects 
at exceedingly high energies, and the approach to the black disc bound in 6-space. Our conclusions are 
presented in Sec. 7. 



2. The GLM Model 

2.1 Single channel model 

The main assumption of the single channel GLM model is that hadrons are the correct degrees of freedom 
at high energy, diagonalizing the scattering matrix. This model [2] fits atot, &ei and B e i well, but fails to 
reproduce the inelastic diffractive final states. This is evident in the relatively better measured SD channel, 
in which the calculated normalization and energy dependence of a s d fail to agree with the experimental 
data [2]. This is not surprising as the input assumption of this class of models, that ^fj is negligibly small, 
is not compatible with the data [10]. 



Experimental data for 
p + p and anti p + p 

R .l= <*„i(s>'<* t „,(s) 



log(s/s ) 



Experimental data for 
p + p and anti p + p 

r d = Ki( s ) + °sd( s ) +°dd( s ) V° tot (s) 



log(s/s ) 



Figure 1: Energy dependence of R e 



Figure 2: Energy dependence of Rd. 



It is instructive to see the ISR-Tevatron experimental data presented in Fig. [I] and Fig. |2[ Fig. [T] shows 

with s , which is compatible with our parametrization (soft 



the power like increase of the ratio R e 
Pomeron like) , and can be reproduced in a single channel eikonal model in the ISR-LHC energy range [11]. 
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Denning Udiff = a s d + Odd-, the ratio Rd = ael+CTdlff which is shown in Fig. |2~L behaves approximately 

J J (7 tot ' 

as a constant of about 0.37-0.38. This is incompatible with the assumed exchange of an unscreened soft 
Pomeron, where, R e \ and Rd are expected to have approximately the same energy dependence. In a single 
channel eikonal model adlff is a small parameter, i.e., diffractive scattering is treated as a perturbative 
effect. As such, this ratio tends to a small constant at low energies, and approaches zero in the high energy 
limit. This behavior is in contradiction with the experimental data shown in Fig. ||. 

In the following we show that the deficiencies of the single channel eikonal model are eliminated in a 
more elaborate two channel model in which diffractive, alongside elastic, scatterings are included in the 
rescattering chain. 

2.2 Two channel model 

The GLM two channel model has been described in our previous publications (see Refs. [1, 12-14] and 
references therein). In this formalism, diffractively produced hadrons at a given vertex are considered as 
a single hadronic state described by the wave function <!'£>, which is orthonormal to the wave function 
of the incoming hadron (proton in the case of interest), < ^h\^D >= 0. We introduce two wave functions 
and ^2 that diagonalize the 2x2 interaction matrix T 

A i,k' =< *fcl T l**' >= A i,k k,v S k ,k', (2.1) 
In this representation the observed states are written in the form 

^ = a^i + /5^ 2 , (2.2) 

^ D = -/3^i + a^ 2 , (2.3) 

where, a 2 + f3 2 = 1. 

Using Eq. fl2.1| ) we can rewrite the unitarity equations in the form 

Im Ai tk (s, b) = \A i>k (s, b) | 2 + G&(a, b), (2.4) 

where G™ k is the contribution of all non diffractive inelastic processes, i.e., it is the summed probability 
for these final states to be produced in the scattering of particle % off particle k. The simple solution to 
Eq. ( |2.4D has the same structure as in the single channel formalism, 

A i>k (s, b)=i(l- exp ) , ( 2 .5) 

GfX(s,b) = l-exp(-Q iik (s,b))- (2-6) 

From Eq. (^1]) we deduce the probability that the initial projectiles (i, k) reach the final state interaction 
unchanged, regardless of the initial state rescatterings, is given by Pf k = exp (— Jlj^s, &)). 

Our presentation, thus far, is model independent. Model dependent parameters are introduced so as 
to obtain explicit expressions for the opacities J7j ^(s, b). As stated in the Introduction, we now replace the 
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simple Pomeron J-pole dynamic input, by an updated interpretation, in which the Pomeron is replaced 
by the description of a partonic system, which goes through a process of saturation when approaching the 
soft, large distance limit [3]. This point of view is compatible with eikonal type models which reproduce 
DIS e-p scattering data over the whole Q 2 range, i.e. from low to high photon virtualities [4,5]. 

Even though we now suggest a more elaborate partonic description for the dynamics of interest, we 
still maintain a simple form for the opacities Oj 

n iik {s,b) = u itk {s)T(s,b). (2.7) 

In the above 

n,k{s) = o^, (-) . (2^ 

and 




The energy term Vi f. is a power in s reflecting a long standing observation [15] that the cross section 
of DIS collisions behaves, at least approximately, as a power in energy. The b-profile Tn- (s, b) are as- 
sumed to be Gaussians. Since the t-space transform of a Gaussian in impact parameter is proportional to 

exp ( — — ) , our parametrization provides a very good reproduction of the small \t\ elastic forward cone 



4 

data, covering more than 95% of the elastic scattering events. A general property of a b-Gaussian profile 
is that the small t region transforms to high b, and high t to small b. The structure of our opacities, as 
presented in Eq. (|2.7|), has the advantage that it is general enough to serve as a parametrization of wide 
class of dynamical models. It is applicable both to a conventional soft Pomeron exchange, as well as a hard 
Pomeron process such as 7 + p — > J/*$> +p [16]. 

Even though we now suggest a more elaborate partonic description for the dynamics of interest, we 
still choose for our radius a Regge-like expression 

Rlk (*) = + R lk + 4 «' Hs/so), = Rl, hk + 4a' ln(-) (2.10) 

so 

in which ■ and a' are fitted parameters. Indeed, the fitted values of a' corresponding to Models B(l) 
and B(2), to be discussed in the next Section, differ from the values typical to Regge phenomenology. Note 
that Rl 2 = 0. 

In general, we have to consider four possible re-scattering processes. For the case of p-p (or p-p) the 
two quasi-elastic amplitudes are equal = 02,1 > and we thus have three rescattering amplitudes: elastic, 
SD and DD. These amplitudes are presented in the two channel formalism in the following form [1, 12, 14] 

a e i(s, b) = i{a*A 1:1 + 2a 2 (3 2 A 1)2 + /3i4 2 , 2 }, (2- 11 ) 

a sd (s, b) = iaP{-a 2 A hl + (a 2 - (3 2 )A 1)2 + p 2 A 2 , 2 }, (2.12) 
a dd = ia 2 p 2 {A hl - 2A h2 + A %2 \. (2.13) 
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It should be stressed that in this approach diffraction dissociation appears as an outcome of the elastic 
scattering of and \&2, the correct degrees of freedom of our model. 

The corresponding cross sections are given by 



a tot {s) = 2 J d 2 blma el (s,b) , (2.14) 

a el (s) = J d 2 b\a el (s,b)\ 2 , (2.15) 

<r s d(s) = J d 2 b\a sd (s,b)\ 2 , (2.16) 

<tm{s) = I d 2 b\a dd (s,b)\ 2 . (2.17) 



In a simplified version we considered a two amplitude model, where a dd is assumed to be small enough 
to be neglected implying that ^2,2 = 2^4^ — A± i. In this model [1] only elastic and SD rescatterings 
are included in the eikonal screening correction. In our past publications we referred to the GLM eikonal 
models according to the number of the rescattering channels considered, i.e. elastic [2], elastic+SD [1] and 
elastic+SD+DD [12]. In retrospect, we consider it more appropriate to define these models according to 
the dimensionality of their amplitude base. We, therefore, call the 2x2 configuration a two channel model, 
making the distinction between its two and three amplitude representations. 



3. Fit to the ISR-Tevatron Experimental Data 

We have studied three models, with different parameterization of Oj which were compared with the global 
experimental data base. The models are based on the general formulae given in Eq. fl2,5| ) - Eq. ( |2.10 ), 



with different input assumptions. Note that the fit has, in addition to the contribution in the form of 
Eq. ( |2.7D , also a secondary Regge sector (see Ref. [1,2]). This is necessary as the data base contains many 
experimental points from lower ISR energies. A study of the vacuum component alone, without a Regge 
contribution, is not possible at this time, since the corresponding high energy sector of the available data 



base is too small to constrain the fitted parameters of Eq. (2/7). Unlike the assumed Pomeron exchange in 
the vacuum channel, the existence of the secondary Regge trajectories is well established both theoretically 
and experimentally. As the goal of this paper is to obtain predictions in the LHC and Cosmic Rays energy 
range, where the contribution of the secondary Regge trajectories is negligibly small, the Regge parameters 
are not quoted in this paper, and will be given in a separate publication. Note that at the Tevatron the 
Regge sector contribution is less than 1%. 

3.1 Model A 

Model A is a two amplitude model, which was considered in Ref. [1] in detail. Its main assumption is that 
the double diffraction cross section is small enough to be omitted from the fitted data base. As we saw, 
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this allows us to express ^2,2 in terms of fi^i and Qi t 2, as such, this model breaks Regge factorization. We 
obtain (see Refs. [1,14]) 

»«,(»,» = i (l - exp - 2^ex P (, -exp (-«))) , (3.1) 



a sd (s, 6) = -ia/3exp ^ 1 - exp ^ , (3.2) 



where, Afi = O12 — ^1,1 • Following Ref. [1], we assume both Qn and AS7 to be Gaussians in b. 



°u ( ^ A 



ni ' i(a « 6)= ^w UJ exp rw^y ' (3 - 3) 



An(s,b) = -^ r (-) (--*). (3.4) 

Note that in this two amplitude model R\ = ^-Rq-ii + 4a'in(s/so) is the radius of Af2(s, b). We have also 
studied a two amplitude model in which both ^11 and fii^ are Gaussians in 6. The output obtained in 
this two amplitude model is compatible with Model A. This is not surprising as a\ 3> crj 1 ( see Table 1). 



3.2 Model B 



In the three amplitude model we do not make any assumptions regarding the values of the double diffraction 
cross sections [6] which are contained in our fitted data base. We use the formulae of Eq. (fj|) and Eq. fl2.10D 
to parameterize the three opacities Qi t i, fii,2 and O2.2, which are taken to be Gaussians in b. If we assume 
the soft Pomeron to be a simple fixed J pole the knowledge of two out of those three opacities determines 
the third. We denote this option Model B(l). In this model cr® 2 = \J a i 1 x a 2 2 ■ ^ s we shall see, the fit 
corresponding to Model B(l) is not satisfactory. We have, also, studied Model B(2) in which a factorization 
of a® k is not assumed. Accordingly, <tJ l5 o"° 2 and a 22 are independent fitted parameters of the model. 
This model reflects our approach to the soft interaction as the continuation of the hard processes in the 
saturated soft region. We would like to stress, that this model gives a very good reproduction of the data. 



3.3 Results 

The parameters of Model A are based on a fit to a 55 experimental data points base which includes the 
p-p and p-p total cross sections, integrated elastic cross sections, integrated single diffraction cross sections, 
and the forward slope of the elastic cross section in the ISR-Tevatron energy range. We did not include the 
Cosmic Ray air showers estimated total cross sections in our data base, as they require additional model 
dependent assumptions [17]. As stated, we neglected the reported DD cross section points. The fitted 
parameters of Model A are listed in Table 1 with a corresponding x 2 1 '(d.o.f) of 1.50. 
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Model 


A 





-Ro;l,l 


a' P 


oil 


°2,2 




<2 


A 


0.126 


0.464 


16.34 GeV~ 2 


0.200 GeV~ 2 


12.99 GeV~ 2 


N/A 


145.6 Gel/" 2 


N/A 


B(l) 


0.150 


0.526 


20.80 GeV~ 2 


0.184 GeV~ 2 


4.84 GeV~ 2 


4006.9 GeV~ 2 


N/A 


139.3 GeV" 2 


B(2) 


0.150 


0.776 


20.83 GeV~ 2 


0.173 GeV~ 2 


9.22 GeV~ 2 


3503.5 Ge^" 2 


N/A 


6.5 GeV~ 2 



Table 1: Fitted parameters for Models A, B(l) and B(2). 



The fits to Models B(l) and B(2) are based on an updated data base which includes the data base 
used for Model A, plus 5 double diffraction cross sections data points [6]. In Table 1 we present two sets 
of fitted parameters for Model B. The factorizable Model B(l) does not give a good reproduction of the 
data, with a x 2 /(d-o./.)=2.30. Part, but not all, of this large x 2 is contributed by the DD data, where the 
model predictions are significantly below the experimental points. Model B(2), with a x 2 /(d.o.f.) = 1.25, 
provides a very good reproduction of the data base. Its seemingly high x 2 reflects the poor quality of the 
published SD data points. The model provides a very good reproduction of the DD data points. 




200 400 600 800 1000 1200 1400 1600 1800 

W(GeV) 

Figure 3: UA4-Tevatron energy dependence of a tot 
in Model B(2). 




200 400 600 800 1000 1200 1400 1600 1800 

W(GeV) 

Figure 4: UA4-Tevatron energy dependence of a e i, 
a s d and add in Model B(2). 



Model A and Model B(2) give compatible reproductions of our data base without the DD points. To 
this end, we refer to the figures of Ref. [1] and Fig. |3] and Fig. Qof this paper. The calculations of B s d and p 
were executed with the fitted parameters of Model B(2). Neglecting the Regge contribution, we reproduce 
the higher energy experimental data points obtained by UA4, CDF and E710. 

Obviously, factorization violations in the coupling input of Model A and Model B(2) are very interesting 
and support our partonic saturation approach. We hope that our present results will encourage additional 
studies of the intriguing relations between, what are traditionally called, soft and hard Pomerons. 
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Figure 5: Total, elastic and diffractive cross sections in Model B(2). 
4. Predictions for LHC and Cosmic Ray Energies 

Model B(2) cross section and slopes predictions for LHC and Cosmic Ray energies are summarized in 
Fig. ||, Fig. ^ and Table 2. At W=14 TeV (LHC energy) our predicted cross sections are: a tot = 110.5 mb, 
a e i = 25.3 mb, a s d = 11.6 mb and add = 4.9 mb. These predictions are slightly higher than those obtained [1] 
in Model A. The corresponding forward slopes are: B e \ = 20.5GeV~ 2 , B s d = 15.9 GeV~ 2 and Bdd = 
13.5 Gey -2 . One can see that the predicted cross section of the diffractive channels, as compared with the 
elastic cross section, is relatively large and should be considered in the estimates for the background for 
diffractive Higgs production process. The ratio of cr e i/atot at the LHC is predicted to be less than 0.25. 
This is a signature that at the LHC energy the proton-proton (anti proton) elastic scattering amplitude is 
well below the black disc asymptotic bound. We shall elaborate on this point in Sec. 6. 

Checking the predictions presented in Table 2, we note a systematic behavior of the listed cross sections 
from which we learn about the gross features of unitarity, spanning the energies from accelerator range to 
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W(TeV) 



Figure 6: Energy dependence of B e i, B sc i and Bdd in Model B(2). 







(J tot 


Pel 




o~dd 


B e i 




Rd 


"difS 
"el 


TeV 


mb 


mb 


mb 


mb 


mb 


GeV~ 2 








1.8 


73.0 


78.0 


16.3 


9.6 


3.8 


16.8 


0.21 


0.38 


0.83 


14 


101.7 


110.5 


25.3 


11.6 


4.9 


20.5 


0.23 


0.38 


0.65 


30 


115.0 


124.8 


29.7 


12.2 


5.3 


22.0 


0.24 


0.38 


0.59 


60 


128.6 


139.0 


34.3 


12.7 


5.7 


23.4 


0.25 


0.38 


0.54 


120 


143.9 


154.0 


39.6 


13.2 


6.1 


24.9 


0.26 


0.38 


0.49 


250 


162.0 


172.0 


45.9 


13.6 


6.6 


26.5 


0.27 


0.38 


0.44 


500 


181.2 


190.0 


52.7 


14.0 


7.0 


28.1 


0.28 


0.39 


0.40 


1000 


202.7 


209.0 


60.2 


14.3 


7.4 


29.8 


0.29 


0.39 


0.10 


10 11 


3970.0 


1070.0 


451.2 


21.6 


19.5 


109.9 


0.42 


0.46 


0.09 


1.22 10 19 


26400.0 


1970.0 


871.4 


25.5 


27.7 


202.6 


0.44 


0.47 


0.06 


(Planck) 





















Table 2: Cross sections and elastic slope in Model B(2). er£Jf is presented for comparison. 



the Planck mass. 

1) The difference between DL non screened total cross section predictions [15] and GLM is small up to 
the GZK ankle, which is a practical upper bound for Cosmic Ray energies with which we can obtain 
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information relevant to our study. 

2) The predicted ratio of R e \ and Rd are well below 0.5 up to GZK and above. Close to the Planck mass 
the ratios approach 0.5. 

3) From Table 2 we see that while R e \ grows very slowly with energy, Rd is essentially a constant, slightly 
less than 0.4, all through the Cosmic Ray energy range. Accordingly, the predictions for the diffractive 
channels, though increasingly suppressed with energy relative to the elastic channel, cannot be ignored. 

4) Well above the GZK cutoff the ratio R e \ grows more rapidly with energy while (J,iiff/^ei which is less 
than 0.1, diminishes slowly. At the Planck mass we see that -R e /=0.44 and -Rd=0.47. 

5) The above predictions suggest a very slow onset of unitarity constraints on the total and integrated 
cross sections with growing energy. We wish to remind the reader that in eikonal models the asymptotic 
behavior of <Jtot, with a logarithmic accuracy, is ln 2 (s/so), that of a e i is \ln 2 {s / 'sq), while that of (Jdiff is 
only InisjsQ). This issue will be discussed in detail in Sec. 6. 



5. Survival probabilities in the GLM model 

In the following we shall limit our discussion to the survival probability of Higgs production in an exclusive 
central diffractive process, calculated in a three amplitude model. A general review of survival probability 
calculations can be found in Ref. [14]. Our one and two amplitude model calculations have been published 
[12,18]. 

In our model we assume an input Gaussian 6-dependence for both the elastic opacities, specified in 
Eq. (|2.7| ) - Eq. ( 2.10| ), and the similar structured hard diffractive amplitude of interest. The hard diffractive 



non screened cross sections in the (i,k) channel are calculated using the multi particle optical theorem [19]. 
As stated, they are written in the same form as the soft amplitudes 

«5 = "5(*)r&(&), (5.i) 

where, 



rf fe (6) = —L- e ^ . (5.3) 



The hard radii R^ are constants derived from HERA J/ty elastic and inelastic photo and DIS production 
[5,16]. See, also, Ref. [14]. 

The general formulae for the calculation of the survival probability for diffractive Higgs boson produc- 
tion have been discussed in Refs. [12, 14,20]. The structure of the survival probability expression is shown 
in Fig. [?]. Accordingly, 

<|S„| 2 > = ^, (5.4) 
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Figure 7: Survival probability for exclusive central diffractive production of the Higgs boson. 

where, 

N(s) =Jd 2 h d 2 b 2 [A H (s, h) A H {s, 6a)(l - A s «s, (bi + b 2 )))] 2 , (5.5) 
D(s) = Jd 2 b 1 d 2 b 2 [A H (s, h) A H (s, b 2 )] 2 . (5.6) 



A s denotes the soft strong interaction amplitude given by Eq. (|2.5|). The form of Ah(s, b) has been discussed 
in detail in Refs. [12,14]. 

Using Eq. ( 2.11 )-Eq. ( |2.13| ), the integrands of Eq. ( pT5| ) and Eq. (5J3) are reduced by eliminating 



common s-dependent expressions. 

N(s) = J d 2 b 1 d 2 b 2 [AH(s,b 1 )A H ( S ,b 2 )(l-A s (b = bi + b 2 ))] 2 

= J d%d 2 b 2 [(1 - a el (s, b))A%(bx)A%(b2) ~ a sd (s, b) (^(^Af (6 2 ) + A^{bM^{b 2 ) 
-a dd (s,b)A p ^(b l )A^(b 2 )} 2 , (5.7) 

D = J d 2 hd 2 b 2 [^(61)^(62)] 2 . (5.8) 

Following Refs. [14, 20] we introduce two hard 6-profiles 



^) = ^exp^j, (5.9) 



The values B^=3.6 GeV 2 and B^=l GeV 2 have been taken from the experimental HERA data on 
J/$> production in HERA (see Refs. [5,16]). 
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A 



Survival Probability 




Figure 8: Energy dependence of centrally produced Higgs survival probability calculated in Models A, B(l), B(2). 



The calculated survival probabilities have been instrumental in the theoretical interpretation of hard 
LRG di-jets produced at the Tevatron. For details see Ref. [14]. Using Eq. (p^) -Eq. (|5.8|) we calculate 
the survival probability Sfj for exclusive Higgs production in central diffraction. Our results are plotted in 
Fig. g. 

In the following we focus on our LHC predictions based on the above. Sjj for exclusive Higgs production 
in central diffraction has been calculated [14] in the two amplitude Model A. The resulting Sfj = 0.027 
is essentially the same as the predictions of the Durham group [21] and Frankfurt et al. [9]. The main 
conclusion of our present results, obtained in the three amplitude B(l) and B(2) Models, is that opening 
more screening rescattering channels results in a reduction of the calculated value of Sj^. Its LHC value in 
Model B(l) is 0.02 and in Model B(2) it is 0.007. This last result has also been obtained in Ref. [9]. Some 
clarifications concerning the reliability of our Sjj estimates are in order. Our model provides am excellent 
reproduction of the elastic amplitude for b > 1.5 — 2.0/m. However, since we do not reproduce da e i/dt 
outside the forward cone well [22], it would appear that our reproduction of a e i at small b is deficient. 
However, we note that a very good reproduction of da e i/dt at \t\ > 0.2 GeV 2 is obtained in eikonal models 
using dipoles or multi-poles in i-space for the profile function. The difference between these opacities and 
ours is very small at small b values [9, 11]. Accordingly, the error introduced into our calculation of S 2 is 
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estimated to be small. An actual comparison [13] supports this conclusion. 

A realistic estimate of the survival probability is crucial for the experimental discovery of a diffractive 
Higgs at the LHC. An educated guess of its value can be obtained, at an early stage of LHC operation, 
through a measurement of the rate of central hard LRG di-jets production (a GJJG configuration) coupled 
to a study of its expected rate in a non screened pQCD calculation. 

6. The Role of Unitarity at High Energies 

6.1 Unitarity experimental signatures 

As it stands, the experimental support indicating the important role of unitarity considerations is confined 
to the energy range attained by existing accelerators. Experimental signatures associated with unitarity are 
observed in inelastic diffractive channels (soft and hard) at relatively low energies. This is a consequence 
of the unitarity suppression which, at a given impact parameter b, is given by Pi t k(s, b) = e _f2i - fc ( s ' b ). After 
its convolution with the inelastic diffractive amplitude of interest we obtain the survival probability. In 
this context we identify a few significant features: 

1) The experimental data shows a severe moderation of the energy dependence of soft diffraction at ISR 
energies and above, and it increase with energy is much slower than elastic scattering. This observation 
differs from the conventional Regge expectation in which the cross sections of elastic and inelastic diffraction 
are expected to have a similar energy dependence. From Table 2 we see that the ratio <Jdiff/&ei decreases 
monotonically with s from low energies to the Planck mass, where it is negligible. In the ISR-Tevatron 
range this has been observed experimentally [10], and is well reproduced by Model B(2). As we have noted 
our calculations at exceedingly high energy are fully compatible with the asymptotic behavior enforced by 
unitarity. 

2) The predicted small value of Sjj, the survival probability of a given LRG hard inelastic diffractive channel, 
has been verified in the Tevatron studies of LRG hard di-jets in various kinematical configurations. The 
survival probabilities are essential in order to understand both the relative smallness and energy dependence 
of the measured LRG hard di-jets experimental rates, when compared with pQCD estimates. Survival 
probabilities are essential in explaining [21] the large factorization breaking in the rates of LRG hard di- 
jets observed in Fermilab and HERA. The above has been systematically obtained by a few independent 
models, regardless of the method chosen to enforce unitarity. For details see Ref. [14]. 

3) The assessment of Sjj is a fundamental ingredient in the experimental program to discover a Higgs 
boson produced in central diffraction. Its signature is that the diffractive system, observed in the r\ — (j) 
final state lego plot, is accompanied on both sides by a LRG. In this paper we have presented Model B(2) 
predicted value of Sjj = 0.007, smaller than our Model A prediction of Sjj = 0.027. Regardless of our 
particular formulation and estimate, unitarity implies that the experimental rate for a diffractive Higgs, if 
and when it will be discovered, is bound to be approximately two orders of magnitude smaller (order of 
1%) than the pQCD estimate. 

6.2 An amplitude analysis of unitarity at very high energies 
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Figure 9: Impact parameter dependence of Af^ and a e i, a s d, add in Model B(2) at W = 25 GeV. 
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Table 3: Impact parameter be- 
haviour of a e i(s, 6 = 0). 



The basic amplitudes of the GLM model are Ax t \, Ai$, and ^2,2 
whose b structure is specified in Eq. These are the building 

blocks with which we construct a e i, a s d, add, (Eq. (|2.11| )-Eq. ( 2.13|) ). 
The Ai t k amplitudes are bounded by the black disc unitarity bound. 
Checking Table 1, it is evident that ^2,2 is considerably larger than 
tti t i and Oi^. As a consequence, the amplitude ^2,2(3, b) reaches 
the unitarity bound of 1 at very low energies. This blackness ex- 
tends to higher b values with increasing energy. The observation 
that ^2 2(s, 6)=1 at given values of (s, 6) does not imply that the 
physical elastic scattering amplitude has reached the unitarity bound 
at these (s, b) values, see Fig. |9| where we show both the basic am- 
plitudes Ai^, as well as, a e i, a s d and add at W=25 GeV. a e i(s,b) 
reaches the black disc bound when Ai t x(s, b)=Ai^(s, 6)=^42,2(s, b)=l. 
Accordingly, a e i(s,b)=l and a s d(s, 6)=a^(s, 6)=0. This result is in- 
dependent of the fitted value of (3. A fundamental feature of the GLM 
models is that a e i approaches the black disc bound at small b very 



slowly, reaching it at approximately W=30,000 TeV. See Table 3 and Fig. 1C. 

A consequence of Jlj & being central in b, is that Pf k (s,b) is very small at b = and monotonically 
approaches its limiting value of 1, in the high b limit. As a result, given an input (non screened) diffractive 
amplitude which is central in b, the output (screened) diffractive amplitude is peripheral. This important 
feature, as calculated in Model B(2), is shown in Fig. |ll] and Fig. |l2| at the Tevatron and LHC energies. 



- 15 - 




As expected, the diffractive amplitudes have a local minimum at b=0, and a maximum at b ~ 1.4/m for 
the Tevatron, which moves to higher values of b with the growth of energy. This implies a non trivial 
t dependence of dodif /(M^^^/dt in the diffractive channels. These qualitative features are induced by 
Model A, Model B(l) and Model B(2) considered in this paper, even though their detailed behavior, as 
seen in Fig. 13, are not identical. Experimentally, this feature is not easily tested as it requires a fine grid 

The implication of the above is that Ai 2, which describes the wave function scattering of ^2 x ^2) 
has (even at W as low as a few GeV) a black central core in 6-space which expands with energy. However, 
the wave function scattering, corresponding to \&i x 'I'i and *$>i x ^2, are much smaller, well below the 
black disc bound (see Fig. 11). Consequently, the amplitudes A\ y \ and A\p have a different behavior as 



- 16 - 




b in fm b in fm 



Figure 11: Impact parameter dependence of A,^ at the Tevatron and LHC in Model B(2). 



a p] , a , and a HH for Model B(2) 
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Figure 12: Impact parameter dependence of a e i, a s d and add at the Tevatron and LHC in Model B(2). 

functions of (s,b), (see Fig. |lT| ). 

The behavior indicated above at the Tevatron and LHC energies becomes more extreme at ultra high 
energies, when a e i reaches 1 at ever larger b values (see Table 3). When a e i(s,b) = 1, unitarity implies 
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Figure 13: Impact parameter dependence of a s d at the Tevatron and LHC in Models A, B(l) and B(2). 



that at these b values a sc i = add = 0. The diffractive cross sections become, thus, highly peripheral and 
relatively small since they are confined exclusively to the very high b values where the periphery of a e i(s, b) 
is below the black disc bound. This prediction is observed in our calculation at ultra high energies, above 
the GZK ankle cutoff. 



We demonstrate this feature and its consequence at the extreme Planck mass. Fig. |1J clearly shows 
that as the black core of a e \ expands, the difference between Models A, B(l) and B(2), considered in this 
paper, diminishes, being confined to the narrow b tail where a e i(s, b) < 1. Fig. |l|| compares a e i, a s d and add 
in Model B(2). Evidently, the diffractive amplitudes, which are logarithmically suppressed, survive just at 
the high b tail of a e \. The above observations may be of interest in the analysis of Cosmic Ray experiments. 



7. Conclusions 

In this paper we have presented a set of predictions for cross sections, forward slopes and survival probabil- 
ities at ultra high energies focusing on the LHC. There is a relatively small variance between the published 
predictions for cross sections at the LHC. Our SD and DD slope predictions are novel and are important 
for a proper assessment of the soft diffractive background to a diffractively produced Higgs. 

We wish to single out three significant properties of our model: 
1) We do not assume the existence of a soft Pomeron, and we obtain the best fit with a parametrization 
which violates Regge (coupling constant) factorization. This feature is in strong contradiction to the 
elementary features of a soft Pomeron which is defined as a pole in the complex J plane. 
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Figure 14: b dependence, at the Planck mass, of a e i 
in Models A, B(l) and B(2). 
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Figure 15: b dependence of a e i, a s d, add at the 
Planck mass in Model B(2). 



2) The elastic amplitude is much smaller than unity in the accessible energy range. 

3) The survival probability for Higgs diffractive production at the LHC is very small, S 2 H = 0.7%. 

Our results for the cross sections and the slopes at Cosmic Ray energies up to the GZK ankle are 
coupled to our analysis of the role of unitarity at high energies. It may appear that our predictions contain 
two contradicting observations. On the one hand, the centrality in b space of the output elastic amplitude 
forces a peripheral b space behavior on inelastic diffraction [23]. This feature occurs already at low energies, 
and is reflected in the low values of the predicted survival probabilities, and the consequent suppression 
of the soft and the hard inelastic diffraction. On the other hand, we obtain a surprising result that the 
approach of the elastic amplitude to the black disc bound in b space, is so slow that a e i reaches unity at 
b=0 only at approximately W=30,000 TeV in the c.m.s. Hence, the elastic amplitude is below the black 
disc bound all through the measurable Cosmic Ray range. Our result is in sharp contrast to Frankfurt et 
al. [9] who claim that a e \ should reach the black bound already at the LHC. 

An intriguing question is how general our results are; i.e. how reliable is our simple parametrization 
of opacities, and our method of imposing unitarity constraints? The fact that the parameters of our model 
depend on the parameterization of the Regge sector, introduces some, though small, level of uncertainty. 
Another weakness is that the GLM model assumes a Gaussian b dependence of the input amplitudes. This 
deficiency results in our inability to properly predict the diffractive dip structure and its position in da e i/dt. 
As we have seen our model reproduces well [22], the forward cone which contains more than 95% of the 
data. 
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Some additional clarifications concerning our three 
amplitude Model B(2) prediction of S 2 H = 0.7% are m 
place. This value is significantly smaller than our two am- 
plitude Model A prediction of S 2 H = 2.7%. We presume 
that our opacities at small b do not introduce significant 
errors in the estimates of Sjj. This has been explicitly 
shown in Ref. [13] for the simple case of a single channel 
model, where we can assess the b-profile directly from the 
^jf- data. For two channel models the Ai t k amplitudes 
are model dependent and can not be determined directly 
from the data, hence the resulting Sjj is also model de- 
pendent. In spite of this, we do have control over the 
output b-distribution of a e i which can be compared with 
the data, very much like the single channel procedure. 
This is the basis for our intuitive assessment that our 
estimate of Sjj is reasonably reliable, regardless of our 
choice for the b-profile. The results we obtain for Sjj de- 
pend on the values we choose for the elastic and inelastic hard slopes B^ and Bf n . We have determined 
these parameters from the HERA measured [5, 16] J/\l/ photo and DIS production elastic and inelastic 
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Figure 16: The dependence of S 2 at the LHC on 
and B^, the slopes for the hard cross sections. 



slopes. Our sensitivity to these parameters is shown in Fig. |16[ Note that when we change the value of 
Bf n we keep the ratio Vp^jBf n unchanged. Doing so we do not change the cross section of the reaction 
7 + P -» J/V + X(M < 1.6 GeV). 



A possible weak feature of our model is the fact that we possibly do not include Pomeron enhanced 
diagrams, which in the Pomeron calculus are responsible for high mass diffraction. The CDF collaboration 
[24] fitting to a triple Pomeron formalism, with ap(0) = 1 + A, found that at the Tevatron energy: A = 
0.112 ± O.OlO(stat) ± O.Oll(syst). We wish to draw the readers attention to the fact that in the triple 
Pomeron formalism [19], the diffractive mass distribution is divergent if A = 0, and it converges for A > 0. 
Hence, replacing the rich diffractive final states by the one state, we implicitly assumed that the mass 
spectrum of diffraction has a typical finite mass. Based on semi-hard processes large mass diffraction has 
been successfully described in e-p DIS [4,5]. We hope to incorporate these processes in our model as well. 
These semi-hard processes are potentially very important to the calculation of survival probabilities (see 
Refs. [25,26]) and can decrease their value further [26]. 



Regardless of these deficiencies, we believe that a qualitative result similar to ours will be obtained in 
any model in which the effect of unitarity is significant. These are open problems and should be further 
pursued. In particular, since we do not have a theory for non-perturbative QCD, the realibility of our 
calculation should be checked by comparing it with the results obtained for alternate models for soft 
interactions. 
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